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I. INTRODUCTION 
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We point out that in general the Reissner-Nordstrom (RN) charged black holes of general rel- 
ativity are not solutions of the four dimensional quadratic gravitational theories. They are, e.g., 
exact solutions of the R + R 2 quadratic theory but not of a theory where a R a bR ab term is present 
in the gravitational Lagrangian. In the case where such a non linear curvature term is present with 
sufficiently small coupling, we obtain an approximate solution for a charged black hole of charge Q 
and mass M. For Q<tf the validity of this solution extends down to the horizon. This allows us 
to explore the thermodynamic properties of the quadratic charged black hole and we find that, to 
our approximation, its thermodynamics is identical to that of a RN black hole. However our black 
hole's entropy is not equal to the one fourth of the horizon area. Finally we extend our analysis to 
the rotating charged black hole and qualitatively similar results are obtained. 
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In the past few years there has been a considerable interest for black hole (BH) solutions and their properties. One 
of the reasons is the developments in our modern physical theories, in particular in gauge field theories and in string 
theory On the one hand these developments have enriched the palette of known fields that probably play some role 
in Nature and of course couple to gravity because of the universal character of the gravitational interaction. On the 
■ other hand, it is expected that Einstein's theory of gravity gets effectively modified, at some higher energy scale. We 
may have in the gravitational Lagrangian non linear curvature terms and the effective coupling with matter need not 
to be minimal. Also the spacetime dimensionality may be higher than four with the extra dimensions compactified 
by some appropriate mechanism. 

Incorporation of such elements in gravitation turns out to be particularly interesting and has already given rise to 
new and, sometimes, surprising consequences for black hole physics. In particular the discovery of a new family of 
— , , BH solutions for Yang-Mills fields Q has set us free from the so-called "no-hair" conjecture and its related theorems. 
^ We now know that black holes do not necessarily belong to the standard black hole families of solutions, namely the 
y | Schwarzschild, Reissner-Nordstrom (RN) or the Kerr-Newman (KN) ones, and furthermore that in general they are 
5h i not unique. In other words, they may have "hair" . The properties of the Yang-Mills fields that are responsible for the 
J^P| existence of new black hole solutions have been extensively analyzed in Ref. For comments on the way classical 
^ , no-hair theorems can be avoided by the presence of some fields, e.g. the dilaton in effective theories inspired on string 
theory, see Ref. ||. 

With the gravity theories modified one may of course expect that also the thermodynamical properties of black hole 
laws for generalized black holes, may significantly differ from that of standard black holes. However several issues do 
survive. E.g. the first law is still valid for Einstein- Yang-Mills black hole solutions In this paper we are interested 
in a somehow accidental property of standard black holes, sometimes referred in the literature as the area law. It is 
the simple proportionality between the black hole entropy S and its geometric intrinsic quantity, the horizon's area 
A, namely the famous S = A/4, relation (in appropriate units). Interestingly enough, this relationship is still valid 
for some generalized black holes, but there are known cases where this is not true. The underlying pattern, if any, is 
not so clear. For references on cases where this relationship is not valid see 0|. In the same reference an interesting 
entropy formula is derived which, in some sense, systematizes the known results. 

Of particular interest for this paper are cases where the gravitational Lagrangian contains non linear curvature 
terms. For higher dimensionality d > 4 it seems that the area law is generically not valid. It is interesting that for 
Lovelock gravity, the entropy is still given by an intrinsic geometric quantity evaluated at the black hole's horizon ||. 
It differs from the Aj 4 value by a sum of intrinsic curvature invariants integrated over a cross section of the horizon. 
For d = 4 in a generic (Riemann) 2 gravity the area law is valid, but not in a generic (Riemann) 3 theory Q. We 
would like however to point out that the above statements are true in vacuum. As we shall see in this paper, when 
matter couples to gravity, e.g. when an electromagnetic field is present, the area law is not valid in a generic four 
dimensional gravitational theory with (Riemann) 2 terms. 

We consider the 4-dimensional theory described by the Lagrangian 
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C= -^^(R + aR 2 +PR^R^)+£ M , (1) 

l07T 

where R^,R are the Ricci tensor and the scalar curvature with respect to the metric g^ v and g := det\g a {>\. a and f3 
are some coupling constants. We will finally take as the matter part, Cm, the usual electromagnetic Lagrangian. 

It is known that the gravitational part of such a theory can be described in terms of a massless tensorial graviton 
field, a massive scalar field if a ^ 0, and/or a massive tensor field whenever (3 ^ 0, see @, @, @ and references 
therein. From the analysis of the field equations at the linearized level 0], it follows that the "source" of the 
additional massive fields is, for the scalar massive field, the trace T^ M ) of the usual matter stress energy momentum 

tensor T$P ; while for the massive tensorial field the source is some combination of components. Using this fact 
we may understand the motivation and some of the conclusions of this paper: the electromagnetic field, having a 
traceless stress energy momentum tensor cannot excite the scalar massive field of the theory. If, furthermore, (5 — 
then the solution may well not differ from that in Einstein's theory, namely it may be the RN solution. Thus the 
no-hair conjecture seems to hold in this case. In fact a no-hair theorem can be proven JTo| . However, it is immediately 
clear that this cannot be guaranteed if /3 ^ 0. It is the purpose of this paper to study this last case and to show 
that the charged black holes of the theory with Lagrangian given by (Q) do not coincide with the RN or, if rotation 
is included, with the KN metrics. We will have to work in some approximation which however in some limit will 
allow us to study the thermodynamical properties of these charged black hole solutions of the quadratic theory (Q). 
In particular we shall consider only solutions expandable in small alpha and beta. 

The paper is organized as follows. In Section II we present the field equations in the exact and in some approximate 
case. In Section III we proceed mainly with the solution of the approximate field equations in the static and spherically 
symmetric case and discuss the range of the validity of the obtained black hole solution. This allows us to study in 
Section IV its thermodynamical properties. In Sec. V we extend the previous analysis to include rotation. Finally, in 
the appendix we give the proof of the basic relations used for the solution of the field equations. 

We use throughout this paper units in which h = c = G = ks = 1, metric signature ( — h H — h), Riemann tensor 
R-tcd '■— ~®d^ic + ' ' ' an d Ricci tensor R a b := Rg\ c b- Finally Gaussian electromagnetic units are employed. 



II. FIELD EQUATIONS 
A. Exact Field Equations 



The gravitational field equations for the theory in Eq. (|j) read 
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(1 + 2aR)(R^ - -g^R) + ^R 2 g^ 

+ (2a + ^)g^R,p P - (2a + (3)R.^ + f3R^-p' p 

-16- SS M (2) 



-g 5{g^) 



where T M „ is the stress-energy- momentum tensor of the matter fields in the Lagrangian of Eq. (|l|). In this paper we 
will consider as matter field the electromagnetic field F^, coupled minimally to gravity. The field equations for F),„ 
are the Maxwell equations 

y/=gV ll F>»' = d^J—gFn = -4ny/=gf, (3) 

where j u is the electromagnetic current. 

Because of the complexity of these equations it is extremely difficult to obtain exact solutions even in cases of high 
symmetry. Consequently one should be prepared to work with some approximation scheme. We shall now discuss the 
main elements of our approximation method. 



B. Approximate Field Equations 

We concentrate on the case where the quadratic theory is slightly different from Einstein's general relativity. In 
particular let us consider the case where we can neglect the quadratic curvature terms in field equations (||). These 
"linearized" field equations read 
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G a b« 8irT a b 

-(2a + ^)g ab aR + (2a + 0)R, ab - 0OR ab . (4) 

More specifically, this linearization demands that terms like R 2 ,R a bR ab , which are proportional to a, (3 parameters, 
can be neglected on the one hand with respect to R, R a b ones and on the other hand with respect to second derivatives 
of curvature terms which are also proportional to a and (3 parameters. The first requirement, can be achieved by 
choosing sufficiently small values of a and f3. Small values for these parameters is actually what nature has chosen, as 
it is concluded from all the tests general relativity (a = Q = /?) has so far successfully passed. The second requirement, 
however, may not be achievable everywhere in the spacetime since it depends on the particular behavior of curvature 
components. We will discuss it further below, when we will consider the case of charged black holes. 

Assuming that the above requirements are satisfied, the spacetime metric g a b in the quadratic theory in the presence 
of some matter field, is slightly modified from the respective metric g a b in Einstein's theory with the same matter 
field. We can write 



gab = dab + X9ab + Ipab, (5) 

with \\\ <C 1, \^>ab\ "C \g a b\- As it is shown in the appendix, by choosing appropriately Xi^Pab the "linearized" 
gravitational field equations Eqs. (^) are equivalent to the following system of equations 

G^(g a b) = SivT^, (6) 

(°-mo)x = -y2\ (?) 
(□ - m?)W + (l - 5) v m v -W> A a) 

= Wn (t^ - \Tg^\ , (8) 



where 



ml 2 = 6a + 2/3, m^ 2 = -j3. (9) 

Finally should satisfy the condition 

(1>ab ~ ^X9ab)' b = 0. (10) 

In Eqs. @, (Q) and (j|) G^ is the Einstein tensor for the metric g^ v . T^ u is the stress-energy- momentum tensor 
of the matter fields in the original Lagrangian. However in our approximation, where quadratic terms in curvature 
are neglected, the energy tensor needs to be constructed using only the g a b metric. Terms that are omitted in this 
way are (dT a b / dg ^ v )5g ^ v oc 0(curvature 2 ) which although they are linear in \ and ip a b fields they are nevertheless 
quadratic in curvature. For the same reason the metric g^ v in the right hand side of Eq. (|^) can be replaced with 
g^u- Finally, in our appoximation the derivative operators can also be taken with respect to 'g il _ v metric. 

To have positive mass parameters mo, mi we shall hereafter assume the no tachyon constraints 

3a + /3>0, (3<0. (11) 



The above field equations demonstrate the decomposition of the quadratic gravitational theory into a theory with 
a graviton field g a b, a massive scalar field \ and a massive tensorial field ip^v, all of them coupled with the stress 
energy momentum tensor T^ v . In particular one should notice here what we have remarked in the introduction: the 
scalar field \ has as source the trace of 1),„ and thus will not get excited by the traceless T^ v of an electromagnetic 
field. On the other hand, the massive spin-2 field ijj^ in general will. Finally let us remark that there are some issues 
regarding the equivalence of such a decomposition with the initial theory. These will not be discussed in this paper, 
see j8| and references therein. Instead we shall only view these equations as a convenient mathematical reduction step 
for solving in some approximation the initial fourth order derivative theory dealing only with second order derivative 
equations. 
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III. SOLUTIONS 



We are interested in charged black hole solutions of the quadratic theory with the Lagrangian ([!]) . We will consider 
both static and stationary solutions. However, we prefer to discuss in more detail the static case which will allow the 
illustration of our method and the exposition of our main results without particular technical complications. Then, 
in another section, we will extend the analysis to include rotation. 

We consider now the case where the matter field is static and spherically symmetric generated by an electric charge 
Q centered at the origin r = 0. The respective spacetime will be spherically symmetric and static and therefore there 
is a coordinate system where the metric can be written as 

ds 2 = -B(r) dt 2 + A(r) dr 2 + r 2 (d6 2 + sin 2 6d<j) 2 ) , (12) 

Here and 4> are the usual spherical angular coordinates. 

The field equations for the Maxwell field are given by Eq. (||) where the current source has only a time component 
equal to j° = QS(r) / (4irr 2 V AB) . These electromagnetic field equations are easily solved in the background metric 
(O). It turns out that the only non vanishing components of the F^ v tensor are the tr (rt) ones with 

ptr = (13) 



The stress energy momentum tensor for this field configuration is 



T\ = T; = -T e 9 = -T$ = - J^l. (14) 



We now look for solutions to the gravitational equations. 



A. Exact Static Black Hole Solutions 



Case (i): a = 0,/3 = 0. In this case we have Einstein's theory where the solution is known to be unique. It is the 



Reissner-Nordstrom (RN) black hole metric given by Eq. (12) with 



B{r)=A(r)- 1 = l-™ + %. (15) 



r 



where M and Q are respectively the mass and the total charge of the black hole. 

Case (ii): a ^ 0, (3 = Here we can check directly [|2j that the general relativistic RN black hole metric is also an 
exact solution for this theory. The question is whether it is unique. The answer seems to be yes since according to 
Whitt (!(]] a no-hair theorem can be proven for theories with gravitational Lagrangian L = R + aR 2 in vacuum and 
in the presence of electromagnetic matter in the case where the condition R ^ — ^- holds everywhere. 

Case (iii): (3 ^ 0. Now the direct check with the RN solution shows that the RN metric is not a solution, except 
for the vacuum case with zero charge Q = where one solution is known to be the general relativistic Schwarzschild 
black hole. The field equations are very difficult and have not allowed us to find an exact solution. 



B. Approximate Static Black Hole Solutions 



We assume now that there exists some parameter space where one may use the approximate system of second 
order field equations (§|)-(||) instead of the original fourth order derivative equations. After finding solutions to these 
equations one has to check whether they are compatible with the assumptions that we have made when deriving these 
field equations. 

For the stress energy momentum tensor of Eq. ( pT[ ) the solution for g^, in Eq. (|^) is of course the RN metric given 
by Eqs. ( |l2| ) and jlq). Now, in the background of this 'g fll , metric we have to solve Eqs. (0) and (||). Notice, once 
again, that the x field is not excited because T = 0. Taking into account the symmetries of our problem we are left 
with the following set of equations for the field -0^ . 

(V 2 - m 2 )W, w 167rT^, (16) 
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where V 2 is with respect to the RN metric, and the stress energy momentum tensor T^ v is given by Eqs. (Q). 

To proceed we shall consider the case of sufficiently small |/3| parameter such that M m^ 1 (remember m^ 2 := 
— 0). This implies that with r > M we will also have r ^> m7/ . Using this fact we observe that we can obtain the first 
leading term for ^„ very simply. Indeed, in the left hand side of the field equations ( |l6| ) the mass term dominates 
over the derivative terms for r > M ^> m^ 1 . Thus, to a good approximation we can neglect the derivative terms and 
arrive at the following diagonal solution for i/v„ 

2Q 2 2M Q 2 

^tt~ a"! (1 + — , 

2Q 2 2M Q 2 



00 



2Q 2 

-^sm 2 9. (17) 

VflyT 

With a straightforward calculation one can check that the condition (|h]) is satisfied. 

Finally replacing in Eq. (||) the components of the Reissner-Nordstrom metric and of the ip^y field from the 
last equation (|l7]) we obtain the metric for the initial fourth order gravitational theory. This can be written as 



r r 2 J \ m 2 r 4 

f 
m\r l 



-r 2 [l- --fL) (d0 2 + sin 2 0d<p 2 ). (18) 



It is easy to see that in this metric the "Newtonian" gravitational potential is modified. At large radial distances, 
particles will feel an additional repulsive force 4Q 2 / (m 2 r 5 ) besides the force that they feel in the RN spacetime. Note 
however that null geodesies on the 9 = plane coincide with those of the RN metric. 

Up to now we have solved the approximate field equations @-dq) in the case where r > M 3> m^ 1 . However, we 
still have to check whether the conditions for using the approximate field equations are satisfied. With straightforward 
but tedious calculations, we find that these conditions are satisfied for the solution ( |l8[) at sufficiently large distances 
from the black hole's horizon r 3> M. However, if one is interested in a solution that is approximately valid in the 
region from radial infinity down to the horizon, one has still to restrict it to sufficiently small charges \Q\ <C M. 
We will provide here only a dimensional argument which leads to this condition. Remember that for using the 
approximate field equations we needed the curvature squared terms to be negligible with respect to second derivatives 
of curvature terms. For our particular problem the metric is approximately the Reissner-Nordstrom one, and thus a 
typical curvature term will be oc Q 2 /r 4 . Therefore we must have (Q 2 /r 4 ) 2 <C (l/r 2 )(Q 2 /r 4 ). If this condition is to be 
satisfied also near the horizon r « M one gets the restriction Q 2 <§C M 2 . And this is what we assume in the following 
section where we discuss the thermodynamic properties of the approximate solution (|l8|). 



IV. THERMODYNAMIC PROPERTIES OF THE SOLUTION 

As we have seen in the previous section the metric of a charged black hole in the quadratic gravitational theory (Q) 
is not the Reissner-Nordstrom solution. To a good approximation the metric is given by Eq. ( |l8| ) . Its range of validity 
extends down to the vicinity of the horizon provided that M 3> \Q\ and M ^ m^ 1 . This will allow us to calculate 
some of the thermodynamic quantities associated to black holes. In particular, we are interested in the entropy-area 
relationship. Since the metric is known to be spherically symmetric and static, it is easy to employ standard Euclidean 
techniques to obtain the temperature associated to the black hole and an expression for its entropy |L3) . However the 
same results are obtained with a more simple minded thermodynamical treatment [ p"l| ] . We prefer here to follow this 
last method. 

With the argument due to Hawking, which is based on a semiclassical calculation of quantum effects near the 
horizon, one can associate to a black hole a temperature T = k/(2it) where re is the surface gravity of the black hole's 
horizon. Note that for our metric the horizon is at a radial coordinate rn given by 
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r H = M + \/M 2 - Q 2 . 



(19) 



The surface gravity of the horizon ft(rn) is then easily computed 



K(r H ) 



1 dgtt/dr 



y/M 2 - Q 2 



2 V-9tt9rr (M + v/M 2 - Q 2 ) 2 



(20) 



We notice here that k(?"h) has the same value as for the RN black hole of general relativity. Consequently the same 
is true for the temperature. 

This coincidence should not create the wrong impression that everything regarding the thermodynamics-geometry 
relationship for our black hole is the same as for the General Relativity. On the one hand the coincidence itself is 
probably accidental in our approximation and does not survive in the exact solution. On the other hand, even to our 
approximation, the other important geometric quantity of black hole, namely its area A, has not the same expression 
in the RN solution and in our solution given by Eq. (|l^) . In fact the horizon area of our quadratic gravity black hole 
is given by 



A = krtri 1 - 



2Q 2 



(21) 



By replacing into this expression the th of Eq. (|T^) and solving the resulting expression for the mass M of the 
black hole we obtain the fundamental relation of black hole thermodynamics: 



M 



A_ 

871 



A_ 

8tt 



) 2 | 2Q2 



Sir) ~r m 2 



1/2 



(22) 



By differentiation of this equation we obtain the first law of black hole thermodynamics 



dM = 



dM 
~dA 



dA + 



dM 
~dQ 



dQ := T K dS + $ H d<3, 



(23) 



where S and <f> represent, respectively, the entropy and the electric potential at the horizon of the black hole. 

Now the task is the proper identification of the quantities Th, S, $h and Q in (^3|) with the black hole parameters. 
The derivatives appearing in ( p3| ) can be explicitly obtained using Eq. (|2^) . The temperature Th is naturally identified 
with the Hawking temperature T — n(rft) / (2tt) that we have calculated above. Similarly Q is the total electric charge 
of the black hole, i.e. the same quantity that appears in our quadratic black hole metric. On the other hand, the 
electromagnetic potential <&h on the horizon can be computed by direct use of Maxwell equations 



$h = 



)dr = 



''H 



)• 



(24) 



Thus by inserting Eqs. (EOT) and (24) into Eq 



we can, upon integration, identify the black hole entropy S as 



8n 2 Q 2 
m\A 



(25) 



Here we observe that extra terms appear when we consider S as a function of A. Thus in general, the simple 
relationship 5* = ^A no longer holds in quadratic gravitational Lagrangian theories. Notice however that in terms of 
the black hole parameters M and Q, the entropy S has the same value as in RN black hole. The same is true, as we 
have seen above, for Th and an d also for other derived quantities as e.g. the heat capacity since 



r -r ds 



°<3 • 



(26) 



Let us remark once again that this nice coincidence probably happens only in our approximation and will not survive 
in the exact quadratic charged black hole solution. 
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V. EXTENSION TO THE ROTATING CASE 



In this section we will extend our calculations to the case of rotating charged black holes in the quadratic theory 
(|l|). We will work in the approximation where we can use the field equations @-(@)- We will use the same methods 
that we used in the previous sections for the static case. 

The metric can now be written in the form 

ds 2 = g tt dt 2 + 2g tt f,dtd(f> + g rr dr 2 + g^dcj) 2 + geedO 2 . (27) 

As in the static case we can almost immediately solve the approximate field equations. Indeed, g^ v will be the 
known general relativistic Kerr-Newman rotating black hole metric. The scalar field \ wu l n °t get excited since the 
source is the traceless electromagnetic stress energy momentum tensor. Therefore, we are left with the field equations 
for the massive scalar field t/> m „ which are to be solved in the background of the KN metric. Then with the same 
arguments as in static case we find that the solution for i/j m „ is approximately given by 

where T^ N ^ is the energy tensor of the Kerr-Newman black hole. Finally the metric in the quadratic theory is given 
by 

g^-g^ + ^P^- (29) 

Using the metric tensor and the tensor of the Kerr-Newman black hole in Boyer-Lindquist coordinates [[lJJ , 

we obtain the metric g^v 

AJ (+) a 2 sin 2 6 

9tt ~ + ~2 7 (-)' 

2 1 asin 2 ' 3 



a sin 2 I 



.900 = [{r 2 + a 2 ) 2 I(-) - Aa 2 sin 2 0I {+) \ — — , 
P 2 

grr = -r J (+) ' ( 30 ) 



where 



20 2 

J (±) = (l±-£4), (31) 

and p := r 2 + a 2 cos 2 6, A := r 2 — 2Mr + a 2 + Q 2 . Here a is the usual angular momentum per total mass parameter. 

As in the static case we shall take \Q\ <C M to be able to discuss the thermodynamic behavior of the solution. 
Having the metric it is now straightforward to calculate the radius of the event horizon rjj, the surface gravity n(rn) 
and the horizon's angular velocity fin- All of them turn out to have the general relativity values namely 



r H = M + ^AP-Q 2 -a 2 , (32) 



y/M* -a 2 -Q 2 

K(r H ) = o— — ^ , (33) 



= ~2~~~7 — 2 • ( 34 ) 



< H TU, 

Straightforward is also the calculation for the area A of the event horizon 

A = ^\{r 2 l + a 2 ) + -^ 1 ! 



'l'H 



1 + ^-2 -arctg( — ) 



(35) 
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Note that, as in the static case, this is different from the general relativity result. From this equation we can obtain 
M as a function of A, Q and J. 

The first law of black hole thermodynamics now reads 



dM 



dM 
~dA 



dA 



1,J 



T H dS 



dM 

~dQ 
<P H dQ 



A, J 

- n H dj, 



dM 
~dJ 



dJ 



(36) 



where the new quantities not present in the static case are the horizon's angular velocity Qn and black hole's angular 
momentum J. The partial derivatives in ( |36| ) can be explicitly calculated, and the Th, Q, $h are naturally identified 
as in the static case. Hence we are left with a differential expression for the calculation of the entropy S in terms of 
the parameters A, J, and Q. Upon integration of this expression we find 



S(A,J,Q) 



A 
1 



in 2 Q 2 8(4tt) 4 Q 2 J 2 



m{A 



3m 2 A 3 



(37) 



We observe that in general the entropy will not be a function only of the area of black hole, but also of the other 
parameters that characterize it. 

It is interesting to notice that also in the rotational case, the actual value of the entropy, that is the one in terms 
of M, Q, J, is the same as in general relativity since S = Airr^ still holds in our approximation. But probably this is 
not true in the exact solution. 



VI. SUMMARY AND DISCUSSION 



We have looked for solutions representing electrically charged, static and stationary black holes in the theory with 
the quadratic gravitational Lagrangian of Eq. The gravitational field equations which involved fourth order 

derivatives for the metric have been reduced to a set of field equations with second order derivatives for a massless 
tensorial field together with a massive scalar and a massive tensorial field. Solving approximately these equations we 
have seen that the modifications from the general relativistic results come only from the massive tensorial field. The 
scalar field does not contribute since its source is the trace of the matter stress energy momentum tensor which in our 
case is the traceless electromagnetic energy tensor. This is in fact an exact result since, as we have directly checked, 
the Reissner-Nordstrom solution is an exact solution of the initial quadratic gravitational field equations only if the 
massive tensorial field is missing (a ^ 0, /3 = in Eq. (ID)). 

Computing several thermodynamic quantities for our black hole solutions we have found that these quantities retain, 
in our approximation, their general relativistic value. This coincidence will quite probably not survive in the exact 
solution. However it is interesting to note here that for generalized charged black holes in other theories usually one 
finds corrections in the temperature and other thermodynamic quantities, see Ref. [jL5| for a discussion in the context 
of string theory and Ref. |l6| in terms of Kaluza-Klein theories. 

Finally, for our black hole solutions we have found that their entropy does not follow the simple area law S = A/4. 
Such a behavior is known for theories that contain non linear curvature terms in their gravitational Lagrangian and 
have dimensionality higher than four. The respective investigations have been mainly performed with the Lovelock 
Lagrangian 0. To have in this theory non trivial results in four dimensions one should have an effective non trivial 
coupling of the curvature to some other field, see e.g. Ref. in the context of Kaluza-Klein theories. As we have 
shown in this paper, we can have in four dimensions a violation of the area law just by including a minimally coupled 
electromagnetic field in the quadratic gravitational Lagrangian of Eq. ([[]) . 



VII. APPENDIX 



The field equations of the quadratic theory in the approximation where only linear terms in curvature tensors are 
retained, are given by equation ([!]). 

Let us now consider the following metric transformation 

gab — > 9ab = 9ab ~ X9ab ~ ^Pab: (38) 

with x and ipab as scalar and tensorial field perturbation, that is with \\\ -c 1 and ip a b <§C gab- Hereafter, hatted 
quantities will refer to the g ao metric. 
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For the purposes of the paper we will keep only linear terms in \ and V'ab fields. Then the transformation in Eq. 
( |38| ) implies that the Einstein tensor G a b is related to G a b with 

Gab = G ab + ^ D i>ab + 7^(lpp P ;ab ~ i>pa;b' P ~ lppb:a' P ) 

+ \gab{^ m - □ V) 

+ X:ab-{V X )9ab + 0{^\x 2 )- (39) 

Here the derivative operators, being unhatted, are with respect to the g a b metric. 
Substituting Eq. (||) in Eq. ([59]) we obtain (to linear order) 

G a b ~ 8nT a b 

-(2a + ^)g ab aR + (2a + 0)R ;ab - puR ab 

+ ^ D V'ah + ^C0p P ;afc - 1ppa;b'' P ~ ^pb;a' P ) 

+\gab(^ Pq ' pq ~ 

+X;ab-(a X )gab + 0(^ 2 , X 2 )- (40) 



Our aim now is to decrease the order of metric derivatives that appear in the field equations (|4C|). This can be done 
by selecting the tensorial field ip a b and the scalar field x m such a way that the terms on the r.h.s. of Eq. (^) which 
contain higher order than second metric derivatives cancel out. Obviously, with ip a b we can cancel the term □i? a &- 
For this let us choose 

i)ab = 2f3R ab + XRg ab , (41) 

where A is some constant. Then it follows 

Gab = SnTab - (2a + (3 + X)g ab OR +(2a + (3 + X)R, ab 

+ X;ab - 9ab^X- (42) 

In deriving Eq. ( [42| ) we have made use of the Bianchi identity R a b' b — \R;a and while commuting covariant derivatives 
we have neglected quadratic terms in the curvature. 
Now with the obvious choice 

X = -(2a + p + X)R, (43) 

we obtain from Eq. (^) that 

Gab « 87rT ab , (44) 

which is indeed an equation with only second derivatives in the metric g a b- Note that the constant A remains 
unspecified and can be chosen arbitrarily. The field equations for x can be derived from the trace of Eq. ^ 

(6a + 2/3)OR- R = 8ttT. (45) 

With the choice of A 

A = ~, (46) 



and using Eq. (43) we can write Eq. ( |45| ) as 

Sir 

□X - mix = — ^-T, TOq 2 = 6a + 2/3. (47) 
The field tp a b should satisfy a set of equations resulting from the Bianchi identity 



®*-w^w°*y =0 - (48) 
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The field equations for tp ab can be obtained from Eq. (|4|) and Eq. (^l|) 

(□ - m^ ol = l6TrT ab + (A - Aa - (3)URg ab 

+ (^ + l)Rg ab + 2(2a + p)R ]ab . (49) 
P 

With the choice and using Eqs. @), (0) and ([H]) we can write Eq. (|9|) as 

2 

(□ - m?)V> fe " (1 - ^1pp P );ab 

= 16n(T ab - ^Tg ab ). (50) 

Thus our initial field Eqs. (|]) are equivalent to the system of Eqs. (|44|), (0), and ([SO]). In these equations the □ 
operator is with respect to the metric g a b but it can be also taken with respect to g a b in our approximation where we 
keep only linear terms in %, ip a b, (and equivalently, because of Eqs. ( |4l| ) and (|4^), linear terms in R 1 R a b)- This is the 
system of field equations used in the body of the paper. 

Note here that Eq. ( p7| ) is not independent since it is just the trace of Eq. (|5(]). From the physical point of view, 
however, it is interesting since it clearly shows that, at least in our approximations, the massive scalar sector of the 
quadratic theory is missing whenever the trace of the stress energy momentum tensor of matter fields T is zero. 
The case of weak gravitational limit is contained in our results and one may compare with the results of Teyssandier 
which however are only valid in a conveniently chosen coordinate system. 
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